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Abstract 



We introduce the language QML, a functional language for quantum computations on 
finite types. Its design is guided by its categorical semantics: QML programs are interpreted 
ly-^ ■ by morphisms in the category FQC of finite quantum computations, which provides a con- 

' structive semantics of irreversible quantum computations realisable as quantum gates. QML 

, integrates reversible and irreversible quantum computations in one language, using first order 

• strict linear logic to make weakenings explicit. Strict programs are free from decoherence and 

' hence preserve superpositions and entanglement - which is essential for quantum parallelism. 

o\ ' 

o 

1 Introduction 

The discovery of efficient quantum algorithms by Shor |18| and Grover |H] has triggered much 
interest in the field of quantum programming. However, it is still a very hard task to find new 
quantum algorithms. One of the reasons for this situation might be that quantum programs are 
very low level: they are usually represented as quantum circuits, or in some combinator language 
which gives rise to circuits. Here we attempt to remedy this situation by introducing the quantum 
programming language QML, which is based on high-level constructs known from conventional 
functional programming. Though functional (programs are expressions), our language is first order 
and unitary; all datatypes are finite. We will discuss possible extensions in the conclusions, but 
we believe that the approach presented here represents a significant progress towards the goal of 
a natural quantum programming language. 

We present a semantics of our language by interpreting terms as morphisms in the category of 
finite quantum computations FQC, which we introduce here. The FQC semantics gives rise to a 
denotational semantics in terms of superoperators, the accepted domain of irreversible quantum 
computation, and at the same time to a compiler into quantum circuits, an accepted operational 
semantics for quantum programs. 

As an illustration, one of the basic quantum circuits is the Hadamard gate, which is usually 
defined by presenting its matrix: 

had = 7fG -h 

But what does this mean in programming terms? In QML this operation is implemented by the 
following program 

had : Q2 — Q2 
had x = if° x 

then {qfalse | (— 1) qtrue} 
else {qfalse | qtrue} 

We can read had as an operation which, depending on its input qubit x, returns one of two 
superpositions of a qubit. We can also easily calculate that applying had twice gets us back where 
we started by cancelling out amplitudes. 
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An important feature of quantum programming is the possibility to create superpositions which 
have non-local effects. A simple application of this idea is the algorithm in figure ^ to determine 
whether two classical bits, represented as qubits, are equal, which is based on Deutsch's algorithm 
(see |12| . pp.32). It exploits quantum parallelism by querying both inputs at the same time; this 
corresponds to the fact that the expressions if ° a and if° b in our program are not nested. The 
famous algorithms by Shor and Grover rely on a more subtle exploitation of this effect. 

eq : Q2 — ° Q2 — ° Q2 
eq a b — let (x, y) — if°{qfalse | qtrue} 
then (qtrue, if° a 

then ({qfalse | (—1) qtrue}, (qtrue, b)) 
else ({(—1) qfalse | qtrue}, (qfalse, &))) 
else (qfalse, if° b 

then ({(—1) qfalse | qtrue}, (a, qtrue)) 
else ({qfalse | (— 1) qtrue}, (a, qfalse))) 

in had x 

Figure 1: A variant of Deutsch's algorithm 

The reader may have noticed that we do not insist on quantum programs being reversible. 
We will discuss this further in section by comparing classical and quantum computation. It 
turns out that in both cases irreversible computations can be reduced to reversible ones in a 
similar fashion. However, reversibility plays a more central role in quantum computation due 
to the fact that forgetting information leads to decoherence, which destroys entanglement, and 
hence negatively affects quantum parallelism. Thus one of the central features of our language 
is control of decoherence, which is achieved by keeping track of weakening through the use of 
strict linear logic (or just strict logic) and by offering different if-then-else (or, generally, case) 
operators, one that measures the qubit, if, and a second, if°, that doesn't - but which can only 
be used in certain situations. We hasten to add that this intrinsic decoherence is not related to 
the decoherence which is caused by thermal noise in a hypothetical quantum computer. As one 
of the referees has pointed out, control of decoherence is in spirit similar to Reynold's control of 
interference |13j . 

2 Related work 

There are a number of papers on simulating or integrating quantum programming within conven- 
tional functional programming, e.g. Mu and Bird's proposal on modelling quantum programming 
in a functional language Karczmarczuk's use of functional programming to model quantum 
systems [5] and Sabry's proposal to structure embedded quantum programs using virtual val- 
ues [2| . Yet another approach was suggested by Sanders and Zuliani which extends the 
probabilistic guarded command language jl()j by quantum registers and operations on quantum 
registers. 

Peter Selinger's influential paper introduces a single-assignment (essentially functional) 
quantum programming language, which is based on the separation of classical control and quantum 
data. This language combines high-level classical structures with operations on quantum data, 
and has a clear mathematical semantics in the form of superoperators. Quantum data can be 
manipulated by using unitary operators or by measurement, which can affect the classical control 
flow. Recently, Selinger and Valiron jTJ] have presented a functional language based on the classical 
control and quantum data paradigm. 

Selinger and Valiron's approach is in some sense complementary to ours: they use an afflne 
type system (no contraction), while we use a strict system (no weakening). The lack of contraction 
is justified by the no-cloning property of quantum states. However, this does not apply to our 
approach as we model contraction by sharing, not by copying - this is also used in [4]. Indeed, 
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classical programming languages do not implement contraction by copying data, but by sharing 
via pointers. 

Andre van Tonder has proposed a quantum A-calculus incorporating higher order |2(J1 119j 
programs, however he is not considering measurements as part of his language. In |19j he suggests 
a semantics for a strict higher order quantum language based on vector bundles. At the current 
time it is not clear, to us, whether the details of this construction work out. 

Abramsky and Coecke have investigated a categorical semantics for quantum protocols using 
the compact closed structure of the category of finite dimensional Hilbert spaces. They suggest 
that their semantics may be relevant for type systems for quantum programming language. It 
remains to be seen how this relates to our work, since our approach does not exploit compact 
closure. 

All the previous approaches adopt a basically combinatory approach to quantum data: oper- 
ations on quantum data are given by combinators implementing unitary operators. We believe 
that our work is novel in that we are proposing high-level quantum control structures, i.e. we are 
aiming at quantum control and quantum data. 

3 Finite classical and quantum computation 

It is frequently emphasised that quantum computation relies on reversibility because quantum 
physics models reversible processes. This is true, but the same holds for classical computation 
- whether we base our notion of computation on Newtonian physics or Maxwellian electrody- 
namics, the underlying physical processes are reversible for a closed system. Hence we should 
explain irreversible classical computation based on a reversible mechanism. Here, we will develop 
a picture which applies to classical and quantum computation. This makes it easy to identify the 
essential differences and also guides the design of QML which realises structures common to both 
computational paradigms by syntactic constructs established in classical functional programming. 

We introduce the category FQC of finite quantum computations and, for purposes of com- 
parison, the category FCC of finite classical computations 1 . We will interpret QML programs by 
FQC morphisms. It is straightforward to identify a classical sublanguage of QML which can be 
interpreted in FCC; however we will not carry this out in detail. 

Objects of both categories are finite sets, for which we use the letters A, B, C. While classical 
computations are carried out on the elements of those sets, quantum computations take place in 
finite dimensional Hilbert spaces; we write C A for the space generated by A, whose elements are 
functions 2 . A reversible finite computation, that is a closed computational system, is modelled 
by a reversible operation <fi, which is a bijection of finite sets in the classical case, and a unitary 
operator on the Hilbert spaces in the quantum case. We write A — ° U nitary B for the set of 
unitary operators from the space generated by A to the space generated by B, which in the finite- 
dimensional case correspond exactly to norm-preserving linear isomorphisms. The initial state of 
a computation is divided into the input A and the initial heap H, and the final state into the 
output B and garbage G; using cartesian product (x) in the classical and tensor product (<g>) in 
the quantum case. To actually perform a computation we also need a heap initialisation constant 
h, which intuitively sets all memory cells in a defined state, e.g. 0. In the classical case this is 
just an element of the set h S H, while in the quantum case it is an element of the vector space 
h G C H . Such a computational system can be visualised by the following diagram: 





— A 





B — 


h 1 — 1 


— H 




G — 



Note that in the above diagram heap inputs are initialised with a h, and garbage outputs are 

1 FCC may be viewed as a categorical account of a finite version of Bennet's results 0. 

2 C A gives rise to a Kleisli structure, [3], here bind is realised by matrix multiplication. Its Kleisli category is 
the category of finite dimensional vector spaces. 
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terminated with a H. To summarise, given finite sets A, B a morphism (H, h, G, <j>) G FCC A B is 
given by: 

• a finite set of initial heaps H, 

• an initial heap h G H, 

• a finite set of garbage states G, 

• a bijection (f)<EAxHc^BxG, 

while a morphism (H, h, G, cj)) G FQCAf? is given by 

• a finite set H, the basis of the space of initial heaps, 

• a heap initialisation vector h G <C H , 

• a finite set G, the basis of the space of garbage states, 

• a unitary operator <f> G A® H -^unitary B ® G. 

Given two computational systems we can compose them by combining initial and final heaps: 




4>fjoa 

More formally, given the morphisms a and (3: 

a = (H a ,h a ,G a ,(j> a ) e FCC4B 
(3 = {Hp, hp, Gp, $0) G FCC B C 

the composite morphism (3 o a ~ (H, h, G, (f>) is given by: 

H = H a x Hp 

h = (h a ,hp) 

G = G a x Gp 

4> = (G a x 4>p) o (Hp x <j) a ) 

Note that we have omitted some obvious symmetric monoidal isomorphisms for x from the defi- 
nition of <j). We leave it to the reader to construct the identity computation. 
Analogously, given morphisms 

a = (H a , h a , G a , 4> a ) G FQC A B 
p = (Hp, hp, Gp, (j>p) G FQC B C 

the composite (3 o a — (H, h, G, <f>) is given by 

H = H a <x> Hp 

h = h a ®hp 

G = G a ® Gp 

(t> = (G a ®(f)p)o(Hp®<p a ) 

Note that <g> is actually x on the underlying finite sets, since C^igiC 5 ~ C AxB . However, we shall 
use the tensor symbol because we interpret the constructed set as the basis of the tensor product 
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of the associated vector spaces. As in the classical case we omit symmetric monoidal isomorphisms 
for ®. 

We consider two computational systems as extensionally equal if they map the same inputs 
to the same outputs. That is, for FCC, a morphism a = (H, h, G, <fi) £ FCC A B gives rise to a 
function on finite sets Ufcc ot £ A — > B by 



Ax H 



BxG 



A 



B 



Ufcc a 

How do we do this for FQC? There is no sensible projection operation on tensor products. Indeed, 
forgetting a part of a pure state (i.e. a vector of the Hilbert space) leads to a mixed state, which 
is modelled by a density operator 6 € A —o A. This is a positive operator, whose eigenvalues are 
interpreted as the probability that the system is in the corresponding eigenstate. Extensionally, 
quantum computations give rise to completely positive mappings, also called superoperators, see 
□ , pp. 136 or \T5\ for details. Given a = (H, h, G, <j>) S FQC A B we write (j> £ A®H -o supC r B®G 
for the associated superoperator ipp = (f)opo(jy. The heap initialisation vector h £ C H can be 
lifted to a density matrix h £ Dens if by h = \h) (h\. Combining this with the partial trace 



operator tic £ B ® G 



3 su P cr B we obtain U F qc a £ A — ° supe r 
A®H — >- B ® G 



B by 



tTG 



A ■ 



B 



Ufqc a 

in the category of superoperators. 

We say that two computations a, (3 £ FAB are extensionally equal [a = ex t (3), if the induced 
maps are equal; U F a = XJ F (3 where F £ {FCC, FQC}. We define the homsets of FCC, FQC 
as the quotients of the underlying representation by extensional equality. It is straightforward to 
verify that composition respects extensional equality. 

As a consequence of our definition we obtain that the assignment of maps to computations 
gives rise to forgetful functors Ufcc £ FCC — * FinSet and Ufqc £ FQC — * Super. Both 
functors are full 3 and faithful. Hence, our categories FCC and FQC can be viewed just as 
different presentations of FinSet and Super. However, going via FCC and FQC has the benefit 
that we get an implementation of our programs as reversible circuits in the classical case and 
quantum circuits in the quantum case. 

An important class of morphisms arc the ones which do not produce garbage, i.e. where G = 1, 
they give rise subcategories FCC°, FQC° of strict morphisms. All strict maps are isometries, i.e. 
linear maps such that (/ v\fw) = (v\w). However, not all isometries arise from strict computations. 



While FQC and FCC are very similar indeed, the fact that FQC is based on wave mechanics 
enables non-local interaction which is exploited in quantum programming. However, there is also 
a new challenge: the possibility of decoherence. Let 6 £ 2 — > 2 x 2 where 2 = {0, 1} be defined as 
8x = (x, x); which can be easily realised by a CNOT gate. The same implementation gives rise to 
0,2 — ° 0,2 ® Q.2, writing Q2 for the object 2 in FQC. In either case, we can compose this with 
7Ti £ 2 x 2 — > 2 (or 7Ti £ Q2 <8> Q2 — °su P cr Q2) which leads to the following picture: 



h 















1 — 1 


< 


5- 






H 



4>s 



3 In the case of FQC fullness is a consequence of Kraus' decomposition theorem. 
4 This is only due to dimensional reasons, indeed in the domain of our interpretation where all spaces are of a 
size 2™ the functor is full. 
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Clearly, classically we have just denned an inefficient version of the identity m o 5 = I; we copy a 
bit and then throw the copy away. However, the situation is quite different in the quantum case: 
while the implementation is given by the same diagram by replacing classical reversible circuits 
with quantum circuits, the composition is not the identity, it is a measurement operation. That 
is, if we input a pure state like R = {-^ |0) + |1)} the output is a mixed state ^{|0)} + |{|1)} 
corresponding to a random qubit. We have lost the advantage of quantum computation and are 
back in the world of classical probabilistic computations. 

As a consequence of this observation we draw the conclusion that one of the main issues a 
quantum programming language has to address is the control of decoherence. This is somehow 
the opposite of the common view which insists that the no cloning theorem outlaws contraction. 
We observe that the implementation of 6 shares a qubit, but it doesn't clone it; considering R 
again we obtain the EPR state {-^ 1 00) + -j= |H)} after executing only 5. We claim that this is a 
natural explanation of contraction because it is completely uniform in both the classical and the 
quantum case. Indeed, classical functional languages do not implement contraction by copying 
data either. S is strict and therefor maps pure states to pure states. In contrast, operations like 7Ti 
are interpreted by a non-trivial partial trace which introduces decoherence. Hence it is weakening 
which deserves our attention, not contraction. 

4 QML: Rules and semantics 

We introduce here the typing rules and the denotational semantics of QML, the latter gives rise to a 
compilation of QML programs to quantum circuits. The compilation is presented diagramatically, 
implementing it requires some care to make sure that the wires generated by subcomputations 
match as intended. 

4.1 Typing rules 

We will only present the typed syntax of QML, which is based on strict linear logic, the untyped 
syntax is implicit in the typed one. We do allow explicit weakenings annotating a term by by 
a context. This leads to an unambiguous type assignment. Any weakening will be translated 
into the use of a non-trivial partial trace, and hence decoherence in the denotational semantics. 
Another source of decoherence is the use of case, or its special instance if-then-else. We make this 
explicit by introducing two different case-operators: one which observes a qubit and thus leads to 
decoherence; and another which is free of decoherence but requires that we derive that the two 
alternatives live in orthogonal spaces. For this purpose we introduce a judgement t ± u. Another 
novelty of our language is a term-former to create superpositions; we can, for example, write 
{ (qtrue, qtrue) | (qfalse, qfalse) }, to create an EPR state. Note that we are ignoring the factor 

which can be automatically inserted by the compiler. The construction of a superposition also 
requires to show that the participating terms are orthogonal. 

Our basic typing judgements are T h t : a meaning that t has type a under context T. and 
L h° t : a for strict terms. We embed h° in h: 

rh°t:<T 

rht:cr 

To avoid repetition, we also use the schematic judgements r h a t : a where a e { — , o}. We use 
a,r and p to quantify over types, which are generated by 1, a © r, a ® r. Qubits are defined as 
Qa = 1©1. 

T is a context, i.e. a function from a finite set of variables domT into the set of types. We 
write contexts as T — x\ : n, . . . , x n : t„ and use • for the empty context. L, x : r is the context 
L extended by x : r. This operation is only defined if T does not already assign a type to x. 
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For the additive rules, we introduce the operator <g> mapping pairs of contexts to contexts: 

r, x : a® A,x : a = (T <g> A), x : a 

T,x:<j®A = (r <g> A), a; : a if x^ dom A 

• <g> A = A 

This operation is partial - it is only well-defined if the two contexts do not assign different types 
to the same variable. 



4.2 Denotational semantics 

We assign to every type cr the number \<j\ which is the size of a quantum register needed to store 
elements of cr, we also interpret expressions of the form <tUt: 

|1| = 
|(7Ut| = max {\a\, \t\} 
\<j®t\ = |(tUt| + 1 
|,t®t| = \a\ + \t\ 

The interpretation of a type is the FQC object of quantum registers of the right size: [cr] = 
Q 2 H . Contexts T : t„ are interpreted as the tensor product of their components 

[r] = [ti] ® [T2] ® . . . (£> [t„J. A typing derivation T h t : c is interpreted as an FQC morphism 
[i] € FQC [r] [cr], correspondingly, T h° t : a is interpreted as [t] G FQC° [r] [cr]. 

The interpretation of orthogonality is more involved. Given r h° t : a and r' h° u : a where 
|T| = |r'| we interpret a derivation t _L u as a structure (S, /, <?, tp) where S is an object of FQC, 
I G FQC [r] S,g£ FQC [r'] 5 such that [t] = (qtrue® -)o / and [u] = 4>o (qf alse® -) o c/. 

To interpret the operator ® on contexts we define an FQC° morphism Cr a € FQC° [r (g) 

A]([r]®[A]) 



r®A — I I — r 



H r ,. 



— A 



by induction over the definition of T (g) A: If a variable x : a appears in both contexts we have to 
use 5 a E FQC° [cr] ([cr] ® [cr]) which generalises 62, discussed earlier, by applying it in parallel to 
all qubits. All the other cases can be dealt with by applying monoidal isomorphisms. Similarly, 
we define an explicit weakening operator Wr,A € FQC [r ® A] [r]. 

4.3 Structural rules 

We start with the strict variable rule and the non-strict weakening and their interpretations 

rht:cr 

,0 var -3 t weak 

x:crh x:cr r®Aht domA : a 



r®A — 



0W r a 



G t 

Gr-A 



Next, we introduce a let-rule which is also the basic vehicle to define first order programs. 

T h a t : cr 

b ?/, : r 

■let 



A, x : cr h fc u : t 



r ® A \- anh let x = t in u : r 
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orio = o and — otherwise. We leave the condition that A is defined as an implicit precondition 
of this and subsequent rules using <g>. The interpretation of the let-rule is given by the following 
circuit: 

r®A ■ 



H r ,. 



OC 



G t 



Weakenings can affect the meaning of a program. As an example consider: 

y : Q.2 \- let x = y in x^ : Q 2 

This program will be interpreted as the identity circuit, in particular it is decoherence-free. How- 
ever, consider 

y ■ 0,2 \~ let 1 = 1/ in x^ : Q2 

This program is interpreted by a circuit equivalent to the one corresponding to 7Ti 08 shown earlier; 
hence it introduces a measurement. 



4.4 Rules for <g> 

The rules for 1, ® are the standard rules from linear logic. In the case of 1 instead of an explicit 
elimination rule we allow implicit weakening: 



1 - intro r,x:lh a t:a 



1 — weak 



The interpretation of the rules for 1 in terms of circuits is invisible, since 1 doesn't carry any 
information. The interpretation of the rules for <g> is more interesting — the introduction rule 
simply merges the components 



r h a t : a A h a u : r 

T <g> Ah" (t,u) : a®T 



— intro 



r®A — I I — r- 

4>G 



ffr,A h 



G t 
G u 



The interpretation of the elimination rule is similar to the let-rule: 



T h a t : a ® r 
A, x : a,y : t \~ u : p 

r ® A h anfc let (x, y) = t in u : p 



■ ® — elim 



Hr,A h 



P 

G t 



As an example, here is a simple program which swaps two qubits: 

p ■ Q2 ® 0,2 V- let (x,y) =p in (y {} ,x {} ) : Q 2 ® Q 2 
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Again it is important to mark the variables with the empty set of variables. The alternative 
program 

p: Q 2 ®Q 2 V- let (x,y) = pin {y {p \x^) :Q 2 ®Q 2 
would measure the qubits while swapping them. 

4.5 Rules for © 

We represent values in a © r as words of fixed length, as in classical computing. Unfolding our 
type interpretation we have that \<j r] = Q 2 <E> \a U r] where [<r U r] can store a value either of 
\a\ or [r]. To adjust the size we use an easily definable padding operator P ctU t € FQC [erj JctLItJ, 
which simply sets unused bits to 0. 

The introduction rules for © are the usual classical rules for +; note that they preserve strict- 
ness. 

T h Q s : a 

V introi 

T h a inl s : cr©r 




H t -s I- 



r h a i 



r h a inr t : a ( 



intr02 



r - 



C 2 h 



crUr 
Q 2 

G t 

where X is negation. 

We define qtrue^ = inl () x : Q 2 and qf alse x = inr () x : Q 2 . To be able to interpret case 
expressions we introduce a biconditional operation on unitary operators. Given (p,tp e A ^unitary 
B we construct 

°unitary 

Q 2 ®B 

by the following matrix 

[<j)\ip] (true, a) (true, b) = 4>ab 
[<j>\ip] (false, a) (false, b) = ipab 
[4>\ip] (x,a) (y,b) = everywhere else 

As already indicated we have two different elimination rules — we begin with the one which 
measures a qubit, since it is basically the classical rule modulo additivity of contexts. 

T\- c: a®T 
A, x : a h t : p 
A, y : t h u : p 



f® Ah case c of {inl x t \ inr y ^> u} : p 



— elim 



We have ft} e FQC [[A © aj |p] and [«] e FQC [A © r] [/>]. By padding the input we turn 
them into [[ij] , [[zt]] € FQC JA © (er U t)J Jp]. There is no reason why the size of the associated 
heap and garbage should be the same, however, we have that H t + G u = H u + Gt and hence we 
can stretch both maps uniformly to H = H t U H u and G = G t UG u giving rise to </>[[*]] and </>[[«]] 
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of identical dimensions. Hence we can apply the choice operator to construct tp = [4>[itj] ^[[u]]]) 
and with some plumbing we obtain: 



i?r,A 



H t - 



- P 

^ Qi 
^ g 
^ g c 



We can derive if-then-else as 

if b then £ else m = 

case 6 of { inl _ =>■ t \ inr _ u } 
and use this to implement a form of negation: 

mnot : Q 2 — ° Q2 

mnot x = if x then qfalse else qtruc 
However, this program will measure the qubit before negating it. If we want to avoid this we have 
to use the decoherence-free version of case, which relies on the orthogonality judgement: ( 1 ti, 
which is defined for terms in the same type and context r h t, u : A. We will introduce the rules 
for orthogonality later. Intuitively, tlu holds if the outputs t and u are always orthogonal, e.g. 
we will be able to derive qtrue^ _L qfalse^. Hence, we introduce the strict case by: 

T h a c : cr©r 
A, x : a\-° t: p 
A, y : t \-° u : p t _L u 



r <g) A h a case c of 

{inl x t I inr y => u} : p 



- clim° 



It turns out that there is no sensible way to define case if a and r have different sizes. Hence we 
define the orthogonality judgement in a way that it only succeeds, if \a\ = |r| and hence [cr] = [tJ. 

To define the interpretation, we have to exploit the data from the orthogonality judgement 
ft 1 u] = (S,/,.g,V) where V € S ® Q 2 ^unitary [p] and /,<? e FQC° ([A] ® [<r]) 5. We note 
that both morphisms must have the same heap and hence we can construct 

[<f>f\<f> g ] G FQC° (Q 2 ® [A] ® [cr]) (Q 2 ® S). 

Now, the main observation is that we just have to apply the unitary operator <fit± u to make the 
qubit disappear, leading to the following diagram: 



Hf-g I ' ^ ( G c 

Note that we only allow strict terms in the branches of a strict case. In a previous draft of 
this paper we tried to be more liberal, however, this causes problems because the qubit we are 
branching over can be indirectly measured by the garbage. This problem was pointed out by Peter 
Sclinger. 

Using the decoherence-free version if we can implement standard reversible and hence quan- 
tum operations such as qnot: 
qnot : Q 2 — ° Q2 
qnot x = if° x 

then qfalse 
else qtruc 



[<Pf\<i>9 



■s- 

-Q.1— 



— p 
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and the conditional not cnot: 

cnot : Q 2 — ° Q2 — ° Q2 <8 Q2 
cnot c x = if c 

then (qtrue, gnoi x) 

else (qfalse, x) 

and finally the Toffolli operator which is basically a conditional cnot: 
toff : Q 2 -° Q 2 -° Q 2 -° Q2 ® (Q 2 ® Q2) 

iojff c x y = if° c 

then (qtrue, cnoi a; y) 
else (qfalse, (x, yj) 



4.6 Superpositions 

There is a simple syntactic translation we use to reduce the superposition operator to the problem 
of creating an arbitrary 1-qubit state: 

T\-°t,u:o- t±u 
PH 2 + ||A'|| 2 = 1 A,AVO 

r h° {(A)t|(A')«}:a 

= if {(A)qtrue | (A')qf alse} 
then t else u 



The algorithm for the preparation of the one-qubit state to a given degree of precision (which is 
a parameter of the compilation) can be obtained from the one-qubit case of the Kitaev-Solovay 
theorem, see P a g e 616-624. 



4.7 Orthogonality 

Given T h t : a and A h u : a where |A| = |T| we define t _L u by the following rules. The idea of 
t _L u is that there is a boolean observation which tells the two terms apart in every environment. 
The interpretation \t _L u\ = (S, f,g,ip) is defined by induction over the derivations. We present 
here a sound but incomplete formalisation of orthogonality achieving completeness is subject of 
further work. 

r h° t : a r h° u : r 
inl t _L inr u inr t _L inl u 



Here p = a r, we set S = a U r. In both cases / is obtained by interpreting t combined with 
padding and r is given by the interpretation of u and padding. The circuits for ip for these rules 
are given by: 




inl t _L inl u inr t _L inr u 



Let r h° inl t, inl u : a © r and let (£, /, g, t^) be the interpretation of t _L u. From this data we 
are constructing the interpretation of inl t _L inl u as (5, f',g',ip'). We set 5" = 5 © Q2 <S> ^ 
where 7? is the heap needed by inl. We construct /' and g' by applying inl to l,r on the level of 
semantics using the appropriate part of S' as the heap, ip is given by the following diagram: 
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H 

0,1 











=x—>c 









The second rule for inr is done symmetrically. 



t _L u 



(t, v) _L (u, w) (v, t) _L (w, u) 



As above, let T h° (t, v), (u, w) : <j®t and let (S, f, g, tjj) be the interpretation of t _L u to construct 
the interpretation of _L (m, w) as (S' , f , g' , ip') . We set S" = S <£> r and construct /' and 5' 
by pairing with v,w, semantically. 

The definition of ip' is given by the following diagram: 



S 



t 1 !t AqKo = — A^Kl 

{(Ao)t I (Ai)u} i. {(«o)t I (ki)u} 

As before, assume as given the interpretation of t _L u as (5, /, g, V')- We construct the interpreta- 
tion of the conclusion as (S, /, g, ip') where ip' is given as 



Q 2 — 6 — . 



using the rotation (j> e Q2 ^unitary Q2 given by 



Ao Ai 
k ki 



4.8 Programs 

So far we have introduced a language of expressions. It is straightforward to extend this to a 
notion of first order programs. E.g. we consider a program S to be a sequence of function 
definitions of the form FT = t : a, we have to parameterise every judgement by S and require 
that r hs t : a for the definition to be a wellformed extension of X. We also have to introduce a 
rule for function-application which can just be translated into an iterated let-expression. 



5 Conclusions and further work 

We have introduced a language for finite quantum programs which uniformly extends a finitary 
classical language. The classical part of our language may be of interest for its own sake, as it 
introduces a natural way to compile functional terms into space efficient reversible circuits, due 
to no unnecessary garbage. This uniformity is one of the main design principles of our language, 
which, we hope, makes it a natural vehicle to express quantum programming and to develop 
quantum thinking. 

We are currently implementing a compiler for QML in Haskell. The compiler produces a 
representation of quantum circuits which can be simulated (inefficiently, of course) by our own 
simulator or by using a standard simulator for quantum gates. 
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There are other design ideas for quantum programming languages. A potential criticism of our 
approach is that we leave contractions implicit, which is an operation which depends on the choice 
of basis. However, our type assignment system clearly fixes the places where contractions have 
to happen, and moreover, and we believe more importantly, it fixes the places where projections, 
or tracing, is happening. A central feature of any quantum programming language seems to be 
control of decoherence. 

Having noted this, it seems that decoherence is something you always want to minimise. It is 
straightforward to design an inference algorithm which infers weakenings t dom r such that decoher- 
ence is minimised. Maybe this should be the default, which can be overridden if the programmer 
wants to enforce measurement. 

We would like to have an orthogonality judgement which is complete with respect to the 
denotational semantics. One of the referees commented that we would need an inner product 
judgement to achieve this. We plan to explore this proposal in future work. 

The restriction that case is only allowed for balanced coproducts is a direct reaction to the 
comments of the same referee who pointed out that our previous approach, which involved padding 
the data, is problematic. Indeed, this problem seems unfixable; if we branch over Qi £g) Q2 the 
garbage which is created by padding may indirectly measure the qubit we are branching over. 
Consequently, this approach would not be compositional, and hence should be rejected. The 
inability to deal with quantum control over arbitrary coproducts is a consequence of the fact that 
while we deal with quantum data and control, the structure, i.e. the memory allocation, of our 
data is classical. One way to overcome this limitation would be to use an operational semantics 
which employs a quantum memory allocation. Such a semantics would have to exploit an infinite 
state space, and it is questionable whether such a system is physically plausible. Another direction, 
which seems more feasible, would be to index quantum structures by classical values at compile 
time. 

We have some doubts as to whether the understanding of general recursion and partiality in 
quantum programming is essential, because partiality is only interesting for systems with infinite 
state spaces. Moreover, it is not clear how to observe the termination of such a hypthothetical 
quantum system of unknown runtime without disturbing the computation. 

Higher order programming would be a worthwhile addition to reflect the way many quantum 
algorithms are presented: e.g. the Quantum Fourier Transform can be parameterised by a function 
on quantum words. Recently, Selinger investigated this problem |16| and it seems that currently 
no canonical higher order structure on Super is known. We are investigating whether the category 
of presheaves over Super would provide a sound denotational model for higher order quantum 
computation. This semantics would employ Day's construction to interpret tensor products. 

Another line of work is to reap the benefits of the fact that our language uses high level 
constructs, and develop high level reasoning principles for QML programs. To achieve this, our 
next goal is to give a direct translation of QML to superoperators which factors through the FQC 
semantics presented here. This translation will be based on the implementation of superoperators 
using arrows [H] in Haskell 21 j. A direct consequence of this construction is that the translation 
presented here is compositional with respect to the extensional equality. 

In joint work with Sabry and Vizzotto we are currently developing an equational theory for 
QML, an algebra for quantum programming, which is sound and complete, with respect to the 
denotational semantics suggested here. Since the completeness proof relies on inverting evaluation, 
such a proof also gives rise to normalisation; exploiting the approach developed in 3 for a classical 
system. 
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